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O ' Abstract In Noyelles etal. (2008, Astron. Astrophys., 478, 959-970), a resonance 

c/3 , involving the wobble of Titan is hinted at. This paper studies this scenario and its 
consequences. 

', The first step is to build an accurate analytical model that would help to find 

I the likely resonances in the rotation of every synchronous body. In this model, 

(N 

[ I take the orbital eccentricity of the body into account, as well as its variable 

' inclination with respect to Saturn's equator. Then an analytical study using the 



second fundamental model of the resonance is performed to study the resonance 



. of interest. Finally, I study the dissipative consequences of this resonance. 
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I find that this resonance may have increased the wobble of Titan by several 
degrees. For instance, if Titan's polar momentum C is equal to 0.355MRj (M and 
Rj being respectively Titan's mass and radius), the wobble might be forced to 
41 degrees. Thanks to an original formula, I find that the dissipation associated 
with the forced wobble might not be negligible compared to the contribution of 
the eccentricity. I also suspect that, due to the forced wobble. Titan's period of 
rotation may be somewhat underestimated by observers. 

Finally, I use the analytical model presented in this paper to compute the pe- 
riods of the free librations of the four Galilean satellites as well as the Saturnian 
satellite Rhea. For lo and Europa, the results are consistent with previous studies. 
For the other satellites, the periods of the free librations are respectively 186.37 
d, 23.38 y and 30.08 y for Ganymede, 2.44 y, 209.32 y and 356.54 y for Callisto, 
and 51.84 d, 2.60 y and 3.59 y for Rhea. 

Keywords Rotation • Natural satellites • Resonance 
1 Introduction 

As most of the major natural satellites of the Solar System planets. Titan is locked 
in a spin-orbit synchronous resonance, i.e. its rotation period is very near to its 
orbital period around its parent planet Saturn (see for instance Richardson et al. 
2004 ji25:]). This corresponds to an equilibrium state known as a Cassini state. 
Thanks to the Cassini fly-bys. Titan's gravity field is known well enough to study 
the behaviour of Titan's spin around the Cassini state. In particular, we now know 
Titan's oblateness coefficients J2 and C22 (Tortora et al. \29\ ). 
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In a recent paper, Noyelles et al. (|21|. hereafter cited as Paper I) give a first 
theory of Titan's rotation, with 3 degrees of freedom. In that work, an analytical 
approach and a numerical one, more complete, are compared. Moreover, some 
aspects of Titan's rotation are elucidated, especially a likely resonance involving 
the free libration of Titan's wobble. 

In this paper, I first propose an improvement of the analytical model resulting 
in better agreement with the numerical results. Then I study the dynamics of the 
likely resonance, and discuss its implications. 



2 Convergence of the analytical study to the numerical study 

I firstly recall how to obtain the Hamiltonian of the problem, as explained for 
instance in Paper I. 3 references frames are considered: the first one (ei,e2,e3) is 
centered on Titan's mass barycenter and is in translation with the reference frame 
in which the orbital motion of Titan is described. This is a Cartesian coordinate 
system refering to the equatorial plane of Saturn and to the node of this plane with 
the ecliptic at J2000. The second frame (ni,n2,n3) is linked to Titan's angular 
momentum G = Gu}, and the third one (fi , i2,h) is rigidly linked to Titan. 

The first variables that are being used are Andoyer's variables (see Andoyer 
1926 iSl and Deprit 1967 fTOl), which are based on two linked sets of Euler's 
angles. The first set {h,K,g) locates the position of the angular momentum in 
the first frame (ei,e2,e3), while the second one {g,J,l) locates the body frame 
(fi,f2,f3) in the second frame (ni,n2,n3) (see Fig.[T|- 
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plane pehpendiciilar 
to angular momentum 



equatorial plane 



Inertia) plane 



Fig. 1 The Andoyer variables (reproduced from Henrard 1141 ). 

The canonical set of Andoyer's variables consists of the three angular vari- 
ables /, g, h and their conjugated momenta defined by the norm G of the angular 
momentum and two of its projections: 



Unfortunately, these variables present two singularities: when / = (i.e. the an- 
gular momentum is colinear to fj, there is no wobble), I and g are undefined, and 
when K = (i.e. when Titan's principal axis of inertia is perpendicular to its or- 
bital plane), h and g are undefined. That is why we use the modified Andoyer's 
variables: 



I 



L = GcosJ 
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G 



h 



H = Gco&K 



p=l+g+h 




r = ~h 



R = ^ = P{1 - cosK) = 2Psin- 



2 K 
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with Q = G-L = G{1- cos/) = IGsin^ f 

Paper I details how to obtain the Hamihonian of the problem in the modified 
Andoyer's variables. It takes account of the free rotation of the body, and of the 
external torque due to Saturn (see Eq|2|. 



Y\+Yl g2 I 71-/2 „2 

-yj ' \- 



1-/1-72 



Jx + Ji 



(2) 



The coefficients of the Hamiltonian are defined as follows: 



, MRi 



MR]. 



7? = 2C22 — 



52 = -^ 



where /2|^ is Saturn's J2, d the distance Titan-Saturn, do the mean value of d, n* 
the mean motion associated with do. 

X and y are the first two coordinates of Saturn in the reference frame (fi,f2,f3) 
bound to Titan. If we call x,-, y,- and zi the coordinates of Saturn's centre of mass in 
the reference frame (61,62,63) defined above, we have 



X 



y 



-R3{-l)Ri {-J)R3 {-g)Ri {-K)R3 i-h) 



I \ 

X. 



(3) 



with 



^ cos (/> — sin ^ 



7?3((/)) = 



sin (/) cos (/) 
1 



(4) 
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and 



\ 



I 



1 

^i(</')= cos (/) - sin (/) • (5) 
yO sin(/> cos(|) 

At the exact Cassini state, a=/? — A + 7r = 0, p = r + f2 = 0, = and 
T]g = 0, A and ^ being respectively Titan's mean longitude and ascending node in 
the frame (61,62,63). In Paper I, Henrard & Schwanen's model (i.e.. Titan mov- 
ing on a circular orbit with a constant inclination and a constant precession of the 
nodes) has been used to obtain the "obliquity" at the equilibrium K* and the three 
fundamental periods of the free librations around the Cassini state. Then a numer- 
ical integration has been performed over 9000 years with complete ephemeris for 
Titan (TASS1.6, Vienne & Duriez 1995 llBOl ) 



Table 1 Comparison between the analytical and numerical results with = 0.31, adapted 
from Paper I. The differences are given here relative to the numerical values, but relative to the 
analytical ones in Paper I. 





analytical 


numerical 


difference 


K* (rad) 


1.1204859 X 10"^ 


1.25481164 X 10-2 


10.7% 


K* (arcsec) 


2311.1681 


2588.2348 


10.7% 


Tuiy) 


2.094508 


2.09773 


0.15% 


TAy) 


167.36642 


167.49723 


0.08% 


T». (y) 


306.62399 


306.33602 


0.09% 



Tabdl adapted from Tab. 17 in Paper I, gives a comparison between the an- 
alytical and the numerical studies, concerning the mean equilibrium position of 
Titan's angular momentum K* , and the periods of the three free librations u, v and 
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w around the Cassini state. The significant difference in the equilibrium position 
is striking, while the agreement can be considered as good for the periods of the 
free librations. 

I suspect that this difference is mostly due to a too simple analytical model, 
that is the reason why I propose a more complete one. 

2. 1 An improved model 

The problem is the accuracy of the model of Henrard & Schwanen (2004 iflBl ). 
Henrard developed more accurate models in the particular cases of lo (2005 liMI ) 
and Europa (2005 fl5l and ||T61 ). but Henrard & Schwanen's model is still the 
most accurate that could be applied to any synchronous satellite. Such a general 
model will not be sufficiently precise for the purposes of obtaining an ephemeris 
of the rotation of a body, but can be used to check the reliability of a more com- 
plete model, e.g. a numerical model like the one performed in Paper I or SONYR 
(Rambaux & Bois 2004 [23]). 

Table 2 Solution for the variable ze (eccentricity and pericentre of Titan), expressed in terms of 
complex exponentials (from TASS 1.6, Vienne & Duriez 1995 1301 ). 



n 


ampl. (rad) 


phase (°) 


frequency [rad.y ') 


period (y) identification 


1 


0.0289265 


153.988 


0.00893386 


703.30 (^6 


2 


0.0001921 


34.663 


-0.00893386 


703.30 -06 



Tab. |2] and [3] come from TASS 1.6 ephemeris and give the main terms of the 
solutions for ze = eexp(v^— IGJ) (eccentricity and pericentre of Titan) and — 
Fexp(\/— lf2) (inclination and ascending node of Titan) with F = sin j. We can 
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Table 3 Solution for the variable ^6 (inclination and ascending node of Titan), expressed in 
complex exponential (from TASS1.6, Vienne & Duriez 1995 (30)). 

n ampl. (rad) phase (°) frequency (rad.y^') period (y) identification 

1 0.0056024 184.578 0.00000000 oo 0q 

2 0.0027899 355.503 -0.00893124 703.51 06 

3 0.0001312 289.015 -0.00192554 3263.07 *>8 

see that the second term in ze has a very small amplitude compared to the first one. 
On the contrary, the second term in i^g cannot be neglected compared to the first 
term, because its amplitude is half the amplitude of the first term. The amplitude 
of the third term is much smaller. 

For this reason I choose to write Z(, and i^g as follows: 

Z6 = ciexp(^^06) (6) 

= roexp(^<?o) +ri cxp{V^<Pe) (7) 

with ei = 0.0289265, Fq = 0.0056024 and Fi = 0.0027899. The analytical study 
performed in Paper I used ei = and Fi = 0, as in Henrard & Schwanen's work 
(2004 [13]). 

With this model, the way to proceed is the same as in the other studies (Henrard 
2005 |[T4l and ifTSi ). The equations have been developed with Maple software 
to the third degree in inclination and second in eccentricity. The Cassini state 
corresponds to the equilibrium of the Hamiltonian : 
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= — nP + Clfi + n5i{l + 5,)[aisin^^r + a2sinX'cos^rcosp 

+^3 COS 2p (1 - COS 27s:)] + «52 ( 1 + 5s) ( 1 + cos^)^ cos 2a 

(8) 

+b2 sin^(l + cos^) cos(2ct + p) + sin^/s:cos(2a + 2p) 
+b4 sin^:( 1 - cos ^T) cos(2c7 + 3p ) + ^5 ( 1 - cos K)'^ cos (2c7 + 4p )] 
with — 0, p — 0, — and T]^ = 0, the coefficients a,, bi being now: 



a,=-\+3{r^+r,^)-lel (9) 



a2 = 2ro + ^ - ^ron^ + 3roej + - SFo^ (10) 

2io 4 4io 



., = ^^±5^ (1.) 



r;2 5 , 27 , sr^e? 3 



fc3 = ^(ro2+ri2) (14) 



^4 = ro3 + ^rori2 (15) 

^5 does not appear because of the truncation to the third degree in inclination, 
so 1 propose to use the expression given in Paper 1, i.e. : 
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Then, the three fundamental periods of libration around the exact equilibrium 
should be processed the same way as in II13II . 11411 . II15I and Paper I. 

2.2 Results 

I here present the numerical applications of the model, with the initial parameters 
given in Tab|4l 

Tab. [5] gives a comparison between the circular model (originally given by 
Henrard & Schwanen 2004 II13II ). the numerical model from Paper I and the model 
given in this paper, with or without ei, and with or without Fi . 
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Table 4 Physical and dynamical parameters. They are the same as in Paper I, except for ei and 
Fi that have been added in this paper and come from TASS1.6 (Vienne & Duriez 1995 |30|). The 
chosen value of is arbitrary. The assumption of the hydrostatic equilibrium is up to now 
doubtful for Titan, in particular because it would imply J2 = ^€22, and so cannot be used for 
estimating . 1 have chosen this value firstly because it is not physically absurd, and secondly 
because it gives a dynamical system far from the resonance hinted in Paper 1. Later in the paper 
1 change this value, to be able to study the resonance itself. 



n 


143. 924047849 1399 rad.y-^ 


TASS1.6 (30) 


e\ 


0.0289265 


TASS1.6 (30) 


ro 


5.6024 X 10"^ 


TASS1.6 (30l 


n 


2.7899 X 10"^ 


TASS1.6 (30) 




58232 km 


lAU 2000 (26) 




1.6298 X 10-^- 


Pioneer & Voyager (7) 


M 


2.36638 X IQ-'^M^ 


Pioneer & Voyager (7) 


Rt 


2575 km 


lAU 2000 (26) 




3.77747586645 x lO^^.km^.y- 


Pioneer, Voyager + lERS 2003 


h 


(3.15±0.32) X 10"^ 


Cassini (29) 


C22 


(1.1235±0.0061) X 10-5 


Cassini (29) 


c 

MRj- 


0.31 





Table 5 Comparison between the different models, with --^ = 0.3 1 . The relative differences with the results from the numerical model are given in parentheses. 



Model 


K* (arcsec) 


Tu{y) 


Tv (y) 


(>•) 


Circular (131 


2311.1681 (10.7%) 


2.09451 (0.15%) 


167.36642 (0.08%) 


306.62399 (0.09%) 


With ei = 


2454.4627 (5.17%) 


2.09452 (0.15%) 


167.36955 (0.08%) 


306.64198 (0.1%) 


With r, = 


2311.1679 (10.7%) 


2.09670 (0.05%) 


167.38996 (0.06%) 


306.64204 (0.1%) 




2454.4625 (5.17%) 


2.09671 (0.05%) 


167.39269 (0.06%) 


306.64198 (0.1%) 


Numerical [211 


2588.2348 


2.09773 


167.49723 


306.33602 
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We can see that the values of the equilibrium "obliquity" K* and of the first 
fundamental period r„ are significantly improved. In particular, taking Fi into ac- 
count helps to approximate K*, while taking ei into account improves the deter- 
mination of r„. 

2.3 Expressing the free solution 

In this part I use the model described in this paper to explain the main free terms 
obtained numerically in Paper I. 

The main terms of the free solution come from the Hamiltonian : 



its derivation being explained in Paper I. Using the different canonical transfor- 
mations that have been used between the variables {c,p,B,q,P,R,r\q) and 
(m,v,w, C/,y,W), we have: 



(17) 



(18) 



v{t) = Oi^t + VQ,V{t)=Vo 



(19) 



w{t) = oi„t + W(i,W{t) = W(i 



(20) 



then 



xi = ^/2UoU*sm{o)ut + uo),yi 




— cos{(Out + uo) 



(21) 
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X2 = yj2VoV*iin{(Oyt + vo),y2 = W -^cos((»,? + vo) (22) 



X3 = y'lWoW* sin(ft)v,? + wo),y3 = -^/^cos(a)v,? + wo) (23) 



and finally 



sin(t()„? + Mo) -j3y2Vby*sin(a)v? + vo) (24) 



p(r) = av^LW*sin(a)„r + Mo) + (l-aj8)\/2W^sin(a)vr + vo) (25) 



2Wq 

r]q{t) = V cos(cOwf + Wo) (26) 



§^(0 = V2WoW*sin((»v.? + wo) (27) 



P(r) = P* + (1 - ap)J^ cos(a^r + «o) - ay ^ cos(o^r + vo) (28) 



7?(?) = 7?* + jS y ^ cos(a)„r + mo) + y ^ cos(a^r + vq) (29) 

where Uq, Vq and Wq are the constant real amplitudes associated with the 3 fun- 
damental proper modes, and mq, vo and wq are the phases at the time origin (i.e. 
J1980, JD 2444240). The constants a, j8, U*, V* and W* are used in the changes 
of variables and are defined in Paper I, and R* and P* are respectively the equilib- 
rium values of R and P. 
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Comparing Eq. (l24l) -( |29] | to the numerical solutions given in Paper I gives the 
following equations: from P: 



77^=9.9514x10-5 



(30) 



from 7?: 




2.28952572 x 10" 



(31) 



from T]g: 



— = 1.514080315 X 10-^ 
W* 



(32) 



from 



y^lWoW* = 3.10703141 X lO"'* (33) 

from ct: 



^2UoU* = 4.78176461 x 10"^ (34) 



|j3|v'2VoV* = 1.147635 x 10-5 (35) 

and from p : 



(1 - aj3)A/2VbV*= 0.18089837 



(36) 



16 



Table 6 Numerical values of the useful parameters, computed with the model presented in this 
paper. They are dimensionless parameters. 



a 


-7.6905124339 X 10-' 




7.0794875675 x 10"' 


U* 


48.03220934680 


V* 


7062.426811769 


w* 


0.2046085670251 



Table 7 Numerical determinations of Uo, Vq and Wq. 



Equation 


amplitude (dimensionless) 


phase 


[30] 


Uo 


= 2.378323538673 x IQ- 


-7 


uo = 


63.00° 


[Ml 


Uo 


= 2.3802020661 13 X 10" 


-7 


uo = 


63.00° 


[a 


Vo 


= 1.851036650588 X 10" 


6 


vo = 


174.72° 


[35] 


Vo 


= 1.860458412872 X 10" 


6 


I'O = 


174.40° 


m 


Vo 


= 2.316782965743 x 10" 


6 


I'O = 


175.64° 




Wo 


= 2.345263498798 x 10 


-7 


m = 


-51.69° 


m 


Wo 


= 2.359051803911 x 10 


-7 


wo = 


-51.69° 



Using Eq [30l[36l and the numerical parameters given in Tab|6| I computed the 
values of the amplitudes of the free librations Uq, Vq and Wq presented in Paper I 
for -rriTr = 0.31. The results are given in TabjTj 

MRj. ° 

We can see a very good agreement between the different determinations of Uq, 
Vo and Wq, which tends to confirm the agreement between the analytical and the 
numerical study. However, Eq [36] gives a slightly different result from Eq[3T]and 
[35] both for Vb and vq. This might be due to a contribution in p whose period 
is very near T,. and so cannot be separated from v over only 9000 years. We are 
limited to this interval of time because it is the limit of validity of the TASS1.6 
ephemeris, as stated by Vienne & Duriez (1995 [3Qi]). 
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3 The resonance 

As seen in Paper I, Titan's rotation seems to encounter a quasi-resonant state when 
Tyj, i.e. the period of the free libration associated with the wobble /, is near 350 
years. This occurs when -^^^ w 0.35. This part aims at first identifying the reso- 
nant argument, then the associated Hamiltonian, and finally to study the associated 
dynamics. 



3.1 Identification of the resonant argument 

From the quasiperiodic decomposition of t]^ + \/~^^q, a periodic contribution 
whose period is about 351 .7 years is likely to lock the system in a resonance with 
w. Unfortunately, such a period might result from different integer combinations of 
proper modes of TASS1.6, more precisely this might be the period of — 2<^6> 206 or 
06 — ^'6. In TASS 1 .6, the amplitudes associated with 06 and <S'6 are respectively e\ 
and Fi . So, we tried several numerical computations with/without ei , with/without 
Fi , to check when the quasi-resonant behaviour occurs (see Fig|2li. We used -~r — 
0.355, to be closer to the exact resonance. 

These numerical integrations have been performed over 600 years, which is 
sufficient visually to discriminate quasi-resonant behaviour from "normal" be- 
haviour. It appears clearly that the quasi-resonant behaviour is reproduced with 
the model presented in this paper, i.e. with one term in eccentricity and two in 
inclination. So, the argument inducing the quasi-resonant state should be / = 

W + <S>6 - 06- 
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■4.5 -+.'1 -4.3 -4.2 -4.1 -4 -3.9 -4.5 -4.4 -43 -42 -41 -4 -3.9 



(c)ei=o,ri^o (d)e, /o,ri/o 

Fig. 2 Quasi-resonant behaviour of Titan's wobble J, with = 0.355 and different models 
of Titan's orbital motion. 

3.2 The Hamiltonian 

I here intend to express the Hamihonian of the problem, considering the resonant 
argument (or quasi-resonant) as a slow argument that should not be averaged. 
I start from the following Hamiltonian : 

JT = ^ + ^ = ©„[/ + w„V + Wh,W + ^ (37) 

where J(f is the Hamiltonian of the complete system, ^ the Hamiltonian of the 
free solution to the first order expressed in (Eqll7b. and 0^ is the remaining part 
of J(f . So, the Hamiltonian is centered on the exact Cassini state and the canoni- 
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cal variables are {u,v,w,U ,V,W). In order for the resonant argument / to appear 
explicitly, I use the following new set of canonical variables : 

M U 
V V 

e = w+^(,-i^6 = f e = w 

Since this change of variables depends explicitly on time, I should add the 
quantity (^'g — ^(,) & to the new Hamiltonian which then becomes ; 

= o)„U + 0),V + (ft),, + <j>6 - 4) + ^2 (38) 

Considering that / = is the only slow argument, every periodic term (except 
of course the slow one) averages to a constant value and is dropped from the 
Hamiltonian, so that ^ becomes: 

= v/0 + jU0^ + eV20cos0 (39) 

The expressions for i//, jx and e are given in the Appendix. 

This Hamiltonian is the Second Fundamental Model of Resonance of Henrard 
and Lemaitre (1983 lfT2]| ') and we can now study its dynamical consequences. 

3.3 Results 

According to (Henrard and Lemaitre 1983 L12J), the resonance associated with 
the Hamiltonian 1391 ) has two stable equilibria if a critical parameter, 5 = — 1 — 



^, is positive. These equilibria correspond to the roots of the 



4_ 

cubic equation — 3{5 + l}x — 2 — with x — V2R and R 



. Only one 
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of these roots is positive, it is of course the one that is interesting and it corresponds 
to a stable equilibrium. 

Table 8 Values of S and of the forced amplitude of the libration with different values of -f-T. 

MRj 

The values of the other parameters are given in TabH The mean wobble <J> has been com- 
puted using Eg 1401 



c 

MRj. 


5 


Wo (forced) 


<J> 


0.34 


5349.4 


(no real solution) 




0.35 


1909.6 


0.342 


80.368° 


0.355 


189.69 


0.108 


40.702° 


0.3555 


17.705 


0.034 


22.337° 


0.355551 


0.1616 


0.010 


12.034° 


0.35555146967191 


1.6x10""' 


0.009 


11.413° 


0.35555146967192 


-3.3 X 10-^ 


(no resonance) 





Different values of 5 and of the forced amplitude of the wobble are given 
in TabUl We see that, under the resonance, the wobble / is increasing if is 
decreasing. 

The mean wobble < / > is obtained as follows: we have + '7^ = ^ ^in^ j 
from Eq(T] Assuming that G w nC, we get ^q+f]q= 4sin^ |. Moreover, it comes 
from EqEilailthat + ,^2 ^ 2U^ [co?,'^ {(0„t + wq) + W*^- ?.iTp-{(0„t +wq)). Equat- 
ing these two last relations we get 

/ = 2arcsin(^^ ^ \J cos2(ft)„,f + wq) + W*'^ ?,ir? {a)„t + wo)^ Finally < / > is 
obtained in averaging this last relation, i.e. 



< y >= ^ y larcsini^'^ ^Jco&^ix) +W*^sm^{T)^dx (40) 
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Since it is a first-order resonance, it is a strong resonance and the capture into 
it is highly likely if the ratio actually has the values given in TablH which is 
still unknown. 

4 Discussion 

4. 1 Consequences for Titan's internal dissipation 

Here I estimate the effect of a forced wobble on Titan's internal dissipation. I use 
for that Wisdom's work (2004 IITll ). in which the reader can find a more detailed 
explanation, as well as in (Murray and Dermott 1999 [20]). 

Time-dependent tidal distortion of a body leads to internal heating. Let Ut be 
the tide -raising potential. The total dissipated energy is given by integration over 
the whole volume of the satellite (here assumed to be Titan), i.e. 



where p is the density, dV the volume element, and v its velocity. Assuming that 
the body is incompressible (it is a good approximation, see Peale and Cassen 1978 
1221 ') and that the density p is constant, Eql41 [becomes: 



where Ut is the tide-raising gravity gradient potential, Uj is the same potential 
including a time delay A due to the tidal response lags of Titan, hs is Titan's 
displacement Love number, and g the local acceleration of gravity. 
The tide-raising gravity potential is 




(41) 




(42) 
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Ut = ^P2(cosa) (43) 

where P2 is the second Legendre polynomial {P2{x) = - \), a is the angle at 
the center of the satellite between the Saturn-to-Titan line and the point in Titan 
where the potential is being evaluated, R is the distance from the center of Titan 
to the evaluation point, and r is the Saturn-to-Titan distance, 
cos a is given by 

o.s 

cos a = — (44) 
rR 

where o and s are, respectively the orbital position of Saturn and the position of 
the surface element in cartesian coordinates whose origin is the center of Titan. 

The cartesian coordinates of a surface element of the satellite in the body-fixed 
frame (fi,f2,f3) are given by 

So = ^ X (sin cos A , sin sin X , cos 6 ) , (45) 

X being the planetocentric longitude, and d the colatitude. Then s is obtained by 
5 successive rotations, the same as Eq|3] i.e. 

s=R,{-l)Ri{ -J)R3 i-g)Ri (-K)R,{-h)so (46) 

Since the goal is just to obtain the contribution of the amplitude of J, I used K = 0, 
h = 0, I = —wt, and g = {n + w)t. In fact, since K is null, g and h are not defined 
but g + h is. Since the node of Titan does not circulate, it disappears in averaging 
over the periods of the motion. Here a constant value for J, Jo, is considered. In 
fact, / is not constant, so 7o could be the mean value of /, i.e. < J >. 
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With these conventions, o = r(cos/, — sin/, 0), / being the true anomaly. 
For small eccentricity. 



r ^ = a '(1 +ecosnf) 



(47) 



cos/ = cos«f + e(cos2«f — 1) 



(48) 



sin/ = sin«f + esin 2nt 



(49) 



and finally dS is given by 



dS = &mededX. 



(50) 



For /o = and a nonzero eccentricity one obtains (see Wisdom 2004 



dE 
'dt 



2 5 



(51) 



With e = and 70 nonzero, I obtain: 



dE 3 3hGMlR\Jl{n + wf 



dt 2 5 

Assuming = 3hs/5 and replacing A by —1/Qs as in Wisdom [31], I get 



(52) 



dE ^3k, GM'fjl^mn + wy 
dt 2Qs na^ 



(53) 



and finally: 



dE 
'dt 



21 , 3 



^-.^ + -(sm/r + -7o^^--j 



GMinRl 



3./« + W\2l fc, „^^"''t) 



X— /- ^ 



(54) 
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where / > 1 is an enhancement factor that takes a partially molten interior into 
account, and / the obliquity of Titan (very small). 

It is widely assumed that the tidal dissipation inside a synchronous satellite 
depends only on the orbital eccentricity, i.e. the other contributions are assumed 
to be negligible. The ratio K between the contribution of the eccentricity and the 
contribution of the wobble is: 



With e = 0.0289 (i.e. TASS1.6 value [30]), k: = 1 (i.e. the contribution of the 
wobble in the tidal dissipation is as high as the contribution of the eccentricity) 
for Jq = 4.381°). For k: = 0.01, we have Jq = 0.438°, so the contribution of the 
wobble in the dissipation is negligible under this value. In addition to this tidal 
dissipation, the wobble induces Coriolis forces that provoke stress and strains, 
responsible for a non-tidal dissipation. This last contribution is very small and 
difficult to determine, because of its dependence on the internal structure of the 
satellite. The reader can find further explanation for instance in Efroimsky and 
Lazarian (2000 ITT\). 

We note that the expression of the contribution of the wobble in the dissipation 
depends on the frequency « + w instead of n. This can be physically explained by 
the fact that the wobble is bound only to the planet (i.e. is not linked, for instance, 
with Titan's orbital plane), so the wobble added to the spin can appear as a sum of 
two motions, whose frequency is n + w. Since the spin is associated with a period 
of 15.94545 days (TASS1.6 value |30|) and the wobble with aperiod of about 350 
years, the composition of these two motions corresponds to a period of 15.94346 




(55) 
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days. The implication here is that the observed rotation may be a little faster than 
that required for spin-orbit synchronisation. Richardson et al. (2004 |25 1) used the 
period of 15.945421 ± 0.000005 days (according to De Pater and Lissauer 2001 
||9l ) and detected a spin period of 15.9458 ±0.0016 days. With a period of wobble 
of 350 years and this period of spin, the detected period might be 15.94343 days, 
which is consistent with their measurements at the three sigma uncertainty level. 

4.2 Application to the other satellites 

I used the analytical model presented in this paper to compute the equilibrium 
obliquity and the periods of the free librations of the Galilean satellites of Jupiter 
and of the Saturnian satellite Rhea (see Tab. 1 10b, using the dynamical and gravita- 
tional parameters listed in Tab|9] 



Table 9 Dynamical and gravitational parameters associated with the four Galilean satellites of Jupiter and the Saturnian satellite Rhea. The mean motion n, 
eu To, Ti and come from Ll-1 ephemeris (Lainey et al. 2006 (l9]) and TASS 1.6 (Vienne and Duriez 1995 |30|), the radii come from the lAU 2000 
recommendations (Seidelmann et al. 2002 1261 ) that are the same as lAU 2006 (Seidelmann et al. 2007 [27]). Nevertheless, here I have used the value of 
Thomas et al. (2006 l28l ) for Rhea's radius, because it is consistent with the parameters of its gravity field derived by less et. al (2007 1171 ). The references of 
the other parameters are indicated in the Table. For Jupiter, I use Ji = 1.4736 x 10^^ (Campbell and Synnott 1985 18)) and Ri^ = 11492km (Seidelmann et al. 
2002 1261 ). For Saturn, these parameters are given in Tab|4] 





J-1 lo 


J-2 Europa 


J-3 Ganymede 


J-4 Callisto 


S-5 Rhea 




1297.2044714 


646.24512024 


320.76544424 


137.51159676 


508.00931975 


ei 


4.15108 X 10"^ 


9.35891 X IQ-^ 


1.42898 X 10-3 


7.37558 X 10-3 


9.713 X 10-* 


To 


3.14217 X 10-"* 


4.04049 X 10"'' 


1.59327 X 10-3 


3.8423 X 10-3 


2.9705 X 10-3 


n 


9.04178 X 10"^ 


2.20043 X 10-'' 


8.53478 X 10-* 


2.24539 X 10-3 


4.207 X 10-* 




-0.845589 


-0.207903 


-0.0456245 





-0.17546762 


M, 


4.70 X IO^'Mt^ f2j 


2.56 X IQ-^Mq^ iZi 


7.84 X IO-^'^Mt^ 18J 


5.60X lO-^M^^ (U 


4.05841 X IQ-^M^ fTTl 


Rs (km) 


1821.46 


1562.09 


2632.345 


2409.3 


764.4 


Ji 


1.846 X 10-3 |3| 


4.355 X 10-* 


1.268 X 10-* (D 


3.110X 10-^ |4| 


7.947 X 10-* (ni 


C22 


5.537 X 10"'' 13J 


1.315 X lO-'' (U 


3.818 X 10-^ (D 


1.050 X 10-^ 14] 


2.3526 X 10-* fm 


c 


0.377 (3] 


0.347 (U 


0.311 □ 


0.353 111 


0.372 [ni 
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Table 10 Equilibrium obliquity and fundamental frequencies of the free librations for the 
Galilean satellites and Rhea, obtained with the model presented in this paper. The results are 
consistent with the ones given by Henrard for To (2005 fT4l ) and Europa (2005 (Ts)). 



Satellites 


K* (arcsec) 


T 




T 


J-1 lo 


140.07 


13.31 d 


159.16 d 


225.17 d 


J-2 Europa 


1866.29 


52.57 d 


3.59 y 


4.84 y 


J-3 Ganymede 


824.02 


186.37 d 


23.38 y 


30.08 y 


J-4 Callisto 


1720.40 


2.44 y 


209.32 y 


356.54 y 


S-5 Rhea 


1320.98 


51.84 d 


2.60 y 


3.59 y 



TabfTTI gives a comparison with previous studies for lo and Europa. We can 
see a good agreement, except perhaps with the "obliquity" of lo. I think that the 
difference of about 10.8% can be explained by the fact that Lainey's ephemeris 
have been updated between Henrard's study and this paper. In particular, Henrard 
used lo = 3.60 x 10~^ rad as given in Lainey's PhD thesis (2002 ifTSll '). which is 
about 12% higher than the value used in this study. 

These data have been computed assuming that the satellites are in hydrostatic 
equilibrium, especially for the values of -f^^- Since the values of the periods 
of libration depend on jf-^, observing these librations should help estimate this 
parameter and test the hypothesis of hydrostatic equilibrium. 

Unfortunately, the free librations should have been damped by dissipative ef- 
fects, unless the amplitude associated is being forced by a resonance (as could be 
the case for Titan). But no resonance appears clearly in TabI 101 (perhaps except 
between Ganymede's wobble W and the node of Europa, with a regression period 
of 30.22 years). 
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Table 11 Comparison of the periods of tlie free librations of lo and Europa given by different 
models. The results labelled "Henrard" come from (Henrard 2005 1141 ) for lo and (Henrard 2005 
(TsJ) for Europa, the ones labelled "SONYR" come from (Rambaux & Henrard 2005 (24)), and 
the ones labelled "circular model" come from Paper I. 





Henrard 


SONYR 


circular model 


this paper 


lo 










K* (arcsec) 


157 






140 


u 


13.25 days 


13.18 days 


13.31 days 


13.31 days 


V 


159.39 days 


157.66 days 


160.20 days 


159.16 days 


w 


229.85 days 




228.53 days 


225.47 days 


Europa 










K* (arcsec) 


1864.62 






1866.29 


u 


52.70 days 


55.39 days 


52.98 days 


52.57 days 


V 


3.60 years 


4.01 years 


3.65 years 


3.59 years 


w 


4.84 years 




4.86 years 


4.84 years 



5 Conclusion 

The general analytical model given in this paper permits a first 3-dimensional 
description of the rotation of every synchronous body. More particularly, it gives 
the equilibrium position, the three fundamental frequencies of the free librations, 
and the main terms of the free solution. It is applied successfully to Titan, lo and 
Europa, in the sense that the results are consistent with the previous studies. 

Moreover, this paper presents a study of a possible resonance involving Titan's 
wobble. With it, I show that if the capture into resonance occured, then the wobble 
could currently have a forced amplitude of several degrees, and so induce a signif- 
icant internal tidal dissipation. Moreover, the forcing of the wobble could falsify 
the space-based detection of Titan's period of rotation to appear a little faster than 
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it actually is, in particular because Titan's rotation axis would be significantly dif- 
ferent from its pole axis of figure. Cassini results on Titan's rotation should give 
some information about that. 

Even if the system is not in resonance, a quasi-resonant state may have a de- 
tectable effect on the free libration w. 
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A Analytical expression of the resonant Hamiltonian 

I here detail the Hamiltonian given by (EqI39K i.e. 




The coefficients \jf, /x and e are respectively: 



(56) 
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64P*'*W*2 V ^ 
+96roV2w*^ - 228i?*P*ro^ + 48W*Vp* - 21W*'*i?*2 
+96W*'^r^P*^ + 16SW*'^r^R*^ - 34m*'^r^R*P* + 24/?* P* 
- 192/f*I^2p*\y*^ - 15R*'^ - 24iy**P*^ + 48ro2p*^ + 24R*P*W*'^ - 8?*%*^ 
+ (2^ - (^)^)(-96roW*V2 + 84/?*roW*V 
-48roP*2 + 48roP*/?* W*^ + eOR'PTo - 48roP*2w*2) j 
+52 ( - 120/?*^ro2 + 2J?*^W*^ - 16ro^P*^W*^ - 16rQ^R*^W*^ 
+22SR*P*r^ + 4&R*P*W*'* - 21R*^W*'^ + 96P*^W*^r^ + mR*^W*^r^ - 348R*P*W*'^r^ 
-24R*P* +40R*P*r^W*^ + 15R*'^ - 24W*'^P*'^ - 48I^^P*^ - 8J?*P*W*^ + SP*^W*^ 

+(2^ - {^f) ( - 96roW*V*2 + 84/?*roiy*V 
+48roP*2 - srTqP'w^ - eoR'PTo + leroP'^w*^) ) ^ 

(57) 

and 



e = """'^^ , (1 + 4) (5i (8P*2 _ 46R*P*W* 

4VW*P'i ^ 

+4P*2 + 16P*^W* - 14R*P* +20R*^W* + 
+(P*V^^(2F^^)(2iro + 69W* -60W*ro^ - 12ro^)) 

+52(-8P*^-46J?*P*W* 
-4R*^ + 16P*^W* + 14J?*P* + 20R*^W*+ 
+ {P* ^yR*{2P*-R*)) ( - 2 IFo + 69roW* - 60W*I^ ^ + 1 2^0 ^ ) ) ) 



We can see that e contains e\ and Fi in its prefactor, wiiicii confimis the fact that there is no 
resonance when one of the two dynamical parameters is null. 
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